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ABSTRACT 


A deployment is the movement of armed forces from their 
home bases to their strategic locations. The movement of 
these forces usually involves the transportation of 
military personnel as well as equipment and supplies. Ina 
crisis situation, it is essential that the deployment is 
carried out in a expeditious manner. 

This study considers the problem of constructing a 
deployment plan for sealift assets which transport military 
personnel, equipment, and supplies to their designated 
locations in the least amount of time. In the construction 
of such a plan, feasible transportation schedules for each 
asset must be specified. When the number of movement 


requirements is large, the problem of arranging schedules 


for the assets is nontrivial. This thesis, therefore, 
describes an algorithm to generate these schedules. Based 
on several examples, this algorithm is shown to be 


effective and can be used in conjunction with algorithms 


for solving the overall deployment problem as well. 
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I. INTRODUCTION 


A. PROBLEM STATEMENT 

A deployment plan is a collection of schedules, one for 
each of the lift assets (trucks, trains, planes, and ships) 
assigned to the deployment. Each schedule normally 
comprises a list of cargoes arranged in the sequence which 
they are to be picked up from the ports of embarkation 
(POE's) and delivered to ports of debarkation (POD's) by 
the designated asset. 

The generation of schedules becomes a great concern in 
the construction of optimal deployment plans because of the 
enormous number of possible schedules. Even for a 
deployment which consists of only 10 lift assets and 20 
shiploads of cargoes, the number of possible schedules 
exceedS 10 million. Thus, regardless of the type of 
algorithm employed for the overall problem, the ability to 
produce an optimal deployment plan in a timely manner 
greatly depends on the ability to produce a good set of 


feasible schedules efficiently. 


B. BACKGROUND 
In mathematical programming, the idea of generating 
columns for the constraint matrix was first introduced by 


Dantzig and Wolfe [Ref. 1]. Among the first to employ the 


column (schedule) generation (or the Dantzig-Wolfe 
decomposition) technique in the area of deployment are Rao 
and Zionts [Ref. 2] and Appelgren [Ref. 3]. Rao and Zionts 
employ the technique to solve the problem of allocating 
ships to existing schedules, and Applegren applies it to 
solve a ship scheduling problem. 

More recently, Brown, Graves, and Ronen [Ref. 4] and 
Brown, Goodman, and Wood [Ref. 5] use a slightly different 
approach in solving the problem of scheduling oil tankers 
and naval surface combatants. Instead of generating the 
schedules as needed in the manner of the Dantzig-Wolfe 
decomposition, they have to generate all feasible schedules 
apriori because the schedules must satisfy many rules, 
restrictions, and relationships among the ports, ships, and 
cargoes. However, these rules, restrictions, and 
relationships also reduce the number of acceptable 
schedules to a manageable amount, thereby making the 


generation of schedules feasible. 


C. OBJECTIVE 

Most authors (see, Ronen [Ref. 6], Ronen [Ref. 7], 
Brown, Graves, and Ronen [Ref. 4], Brown, Goodman, and Wood 
[Ref. 5], Collier [Ref. 8], and Lally [Ref. 9]) consider a 
deployment with the minimum cost to be optimal. Instead of 
cost, one can also minimize other cost related quantities 
such as total distance eaveled® or the number of lift 
assets required. These quantities, if necessary, can always 
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be converted into a monetary amount. However, this study 
adopts a different optimality criterion which arises in the 
military deployment planning during the period of conflict 
(crisis). In this situation, the speed at which the armed 
forces, equipment, and supplies are deployed to the area of 
conflict is one of the most important factors governing how 
the conflict is resolved in the end. Therefore, it is 
essential that the plan can deliver all cargoes to their 
destinations in the shortest amount of time. 

The objective of this study is to develop an efficient 
algorithm to generate schedules that can be used in 
conjunction with any column generation scheme for solving 
the crisis deployment problem. One such scheme is presented 
in a related Naval Postgraduate School Masters Thesis by 
It. N. R. Lima [Ref. 10}. 

The remainder of this thesis is divided into 4 
chapters. Chapter 2 gives a summary of the mathematical 
formulation of the deployment planning problem. Chapter 3 
describes the schedule generator. Chapter 4 discusses the 
computer implementation and the results from the 
experimentation with a group medium size problems. Finally, 
Chapter 5 summarizes the study and indicates areas for 


future research. 


II. PROBLEM FORMULATION 


In the deployment problem addressed below, all cargoes 
are assumed to be in full shiploads. This implies that the 
asset assigned to pick up a given cargo must deliver it 
before any other cargo can be picked up. Because of this 
assumption, the cargoes are also referred to as "movement 
requirements". Moreover, it is also assumed that other 
data essential to the problem, such as the distances 
between ports, ships' speeds, ships! initial positions, and 
tables indicating the ship/port and ship/cargo 
compatibility, are also given. Then, the problem of 
constructing a deployment plan with the minimum duration of 


time can be stated as follows: 


Indices 

i - movement requirements (cargoes), where 1 = 1,....,M 
and M is the number of movement See wiaemantse. 

j} - ships, where j = 1,..... ,N and N is the number of 
available lift assets. 

k - feasible ship schedules. As mentioned in Section C of 
Chapter 1, the generation of these schedules is the 


topic of this thesis. 


Data 





Sijk = 


tik ~ 


1 if and only if the kth feasible schedule for the 
jth ship includes the ith movement requirement. 
Thus, S3k 1s an M-dimensional vector representing 
the kth feasible schedule for ship j. A schedule 
for ship j is considered feasible if ship j is 
compatible with the cargoes it must pick up and 
with the ports it must visit. 


the completion time of schedule S4k- 


Decision variable : 


X4k = 1 if and only if the kth schedule for ship j is 
selected for the deployment. 
Problem Pl 


Min { max ( = tik Xjk : J) = 2), Nee 
k 


J 
Subject to 
N 
= = Sisk Xk Zed ff OL dee — sale yM (1) 
j=1 k 
ZB X4k <1, forj =1, ,N (2) 
k 
Xjk = 020% 1. (3) 


In the above formulation constraint (1) forces each 


movement requirement to be picked up by at least one ship, 


while constraint (2) allows each ship to sail at most one 
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schedule. A set of x's satisfying constraints (1), (2), and 
(3), if it exists, would specify a deployment plan to 
Support all the requirements of the deployment using the 
given assets. 

In a related thesis, Lima [Ref. 10] obtains a solution 
to Problem Pl by parametrically solving a simpler problem 
several times. This simpler problem, Problem P2(T) below, 
is called the feasibility-seeking problem and is a variant 
of Problem Pl. In the feasibility-seeking problem, one 
assumes that a maximum time for deployment, 1, is given. 
Then, the task is to find a set of schedules with 
completion time less than or equal to 7 that satisfies the 
constraints on Problem Pl. To state the problem 


mathematically, let 


Auxiliary Variables 


W; = 1 if and only if movement requirement 1 is not 


delivered, and 


Auxiliary Data 


Kj (7) 


the index set of feasible schedules for ship j 


with the completion time less than or equal to 


K5(T) =—Jils Ks S4k is a schedule and ti. = 7 }. 
Then, the feasibility-seeking problem can be stated as 


follows: 


Problem P2(rT) 


Min 2 w; 


a 
subject to 
N 
2 2 Sijk X43 + Wi Zeeebor- ir = 1, ,M (4) 
J=1 keK3 (1) 
>e X43 Sel ete) <=), ,N (5) 
k 


ie 0 “Or aie, (6) 


Associated with each constraint in (4) there is a dual 
variable u; (s0) which is referred to as the “cargo duals" 
in the subsequent chapters. Similarly, associated with each 
constraint in (5) there iS a dual variable v5 (20) which 
will be referred to as the "ship dual". 

Observe that for a given value Tf, if the optimal 
objective function value of P2(T) 1s zero, then the upper 
bound for Problem Pl is Tt. On the other hand, if the 
optimal objective function value of P2(T) is positive, then 
there is at least one undelivered shipload of cargo. 
Therefore, tr must be increased in order to obtain a 
feasible plan with the same number of assets. By varying T 
and resolving Problem P2(T) in a systematic manner, one can 
obtain a solution to Pl. 

Although Ks (T) is a subset of the feasible schedules, 
it is still too large to enumerate. However, one can avoid 
this enumeration of schedules in the set ay, applying 
the Dantzig-Wolfe decomposition technique to Problem P2(rT) 
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which produces a master and a subproblem. The master 
problem has the same form as Problem P2(T) and the 
subproblem generates schedules (or columns) which belong to 
the set K;(T) and have negative reduced costs. 

The following chapter describes an efficient algorithm 
for the subproblem. The details of solution techniques and 


strategies for the master problem are presented in the 


related thesis by Lima [Ref. 10]. 


III. A SCHEDULE GENERATOR 


A. THE SUBPROBLEM 
The subproblem resulting from the application of the 
Dantzig-Wolfe decomposition to Problem P2(1T) can be stated 


as follows: 


Problem S1(j) : 
k(j) = arg min { Vj + 2 Sijx Uy t tox £7} 
k i 

As previously defined in Chapter 2, uj; and v4, are the cargo 
and ship duals of Problem P2(T). For a given deployment 
plan, a negative cargo dual, uj, indicates that the ith 
shipload is undelivered, and its magnitude can be regarded 
as the profit obtained from delivering the cargo to its 
destination. The ship dual, V3, can be interpreted as the 
cost of using ship Jj. Therefore, the problem S1(j) would 
find a new schedule which utilizes ship j in the most 
profitable manner. 

Moreover, Problem S1 is indexed by j to indicate that 
the generated schedule is for ship j. In theory one can 
solve one subproblem for a particular ship j, or solve N 
subproblems, one for each ship. To insure convergence only 
one schedule (column) with a negative reduced cost needs to 
be added to the master problem during each cycle of the 


decomposition process. 


If the set of schedules Sj, is available apriori, 
solving Problem S1(}) is merely a matter of selecewna one 
schedule from the set Kj;(7T). Otherwise, Problem S1(j) 
resembles the vehicle routing problem [Ref. 11] with the 
exception that, 


(1) the ship schedules do not have to. start and 
terminate at the same port, 


(2) the shiploads of cargoes must be delivered before 
any other shiploads can be picked up, and 


(3) the schedule must be completed before time 17. 
These three additional conditions greatly increase the 
complexity of Problem S1()j). 

As stated, Problem S1(j) would always produce a 
schedule with the most negative reduced cost for each ship 
j}. Moreover, considering the complexity of the problem, it 
may not be advantageous to solve Problem S1(j) to 
optimality. Moreover, the convergence of Dantzig-Wolfe 
decomposition only requires that the reduced cost of the 
new schedule (column) be negative. Thus, one only needs to 
search for a schedule, S34,, which satisfies the following 


two conditions: 


V4 = » Sijk Uj < 0O (7) 
aL 


t4k Ss sT (8) 
where condition (7) ensures that the schedule has a 
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negative reduced cost and (8) guarantees that the schedule 
can be completed within the allowable time, T. 

In practice, entering columns which merely satisfy 
conditions (7) and (8) does not generally lead to an 
efficient algorithm. In the early iterations of the 
Dantzig-Wolfe algorithm, the ship duals, vj, tend to be 
zero and the cargo duals are mostly negative. Thus, a 
schedule which picks up only one cargo with a negative dual 
would easily satisfy both conditions (7) and (8). However, 
such a schedule would be considered a bad schedule if the 
ratio of cargoes (movement requirements) to the number of 
ships is much larger than 1. For example, if the ratio is 
four then on the average one would expect each ship to 
complete four movement requirements. So, a schedule that 
picks up only one cargo is unlikely to be included in the 
optimal plan. The algorithm described below will search for 
one or more schedules which pick up a minimum number of 
cargoes as well as satisfying conditions (7) and (8). This 
minimum number of cargoes is calculated based on the ratio 


of cargoes to ships. 


B. AN ALGORITHM WITH THE DEPTH-FIRST SEARCH STRATEGY 

A schedule which satisfies conditions (7) and (8) and, 
at the same time, picks up the minimum number of cargoes is 
considered to be a "good" schedule. The algorithm below 
searches for one or more good schedules in the same manner 
as the branch and bound algorithm, with the depth-first 
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strategy, searches for an optimal solution for an integer 
program. 

Basically, the algorithm tries to generate a sequence 
of numbers which is called a path. With the exception of 
the number zero, which represents the ship's initial 
position, the numbers in the path represent the cargo 
number and the order of the numbers represents the sequence 
in which the cargoes are to be picked up and delivered. 
Thus, a schedule or column for the master problem can 
easily be constructed from this path. 

To illustrate the principle underlying the algorithn, 
consider a deployment problem with 3 movement requirements 
numbered 1 to 3. Figure 3-1 depicts a tree which represents 
all possible schedules for a ship. The node number in the 
tree corresponds to the movement requirement number except 
for the root node, which iS numbered zero and corresponds 
to the ship's initial position. By traversing this tree 
starting at the root node, one would trace out a path from 
which a schedule can be constructed. For example, consider 
path 1 in Figure 3-1. This path corresponds to the 
following sequence of numbers: 0-1-3-2 which in turn 
corresponds to a schedule for the ship to complete 
movement requirements 1,3, and 2 in sequence. Similarly, 
path 2 corresponds to a schedule for the ship to complete 


movement requirements 2 and 3 in sequence. 
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Initially, the algorithm discards all movement 
requirements which are either incompatible with the ship in 
consideration or have a positive dual variable. This 
incompatibility is due to either ship/port or ship/cargo 
incompatibility. Assume that the remaining movement 
requirements are numbered as 4, 6, 7, and 9. Then, starting 
at node zero with the current path containing only the 
number 0, the algorithm creates nodes 4, 6, 7, and 9 and 
attaches them to node zero (see Figure 3-2). From node 
zero, the algorithm considers adding a new movement 
requirement to the current path, say node 4; this is called 
branching. At node 4, the algorithm checks to see if the 
schedule which contains movement requirement 4 can be 
completed in 7+ days. If the answer is yes, the algorithm 
then includes node 4 into the current path, creates nodes 
corresponding to movement requirements which are not 
members of the current path, and attaches them to node 4. 
These nodes are node 6, 7, and 9. At this point, the 
algorithm would branch to, say, node 6 and at node 6 it 
checks to see if the schedule which contains movement 
requirements 4 and 6 can be completed in +r days. If no, the 
node is not added to the current path. The algorithm would 
then check if the current path, 0-4, corresponds to a 
schedule with negative reduced cost and picks up the 


minimum number of shiploads of cargo. If no, node 6 is 
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"fathomed" and the algorithm would branch to “eurer 
unfathomed nodes in the depth-first manner, which in this 
case would be node 7 or node 9 in the second level. If yes, 
then a good schedule is found, and the algorithm would stop 
if only one good schedule is desired. If more schedules are 
required, the algorithm would store this good schedule, 
fathom the node, and branch to other unfathomed nodes to 
obtain the desired number of schedules or until all nodes 
in the tree are fathomed. 

Below, we formally state the algorithm discussed above. 
It employs the abstract data type called STACKS [Ref. 12] 


to dynamically generate the search tree and conserve memory 


requirements. 
Algorithm 
Input : 
J = the index of the ship being considered. 
T = the maximum allowable completion time for any 
schedule. 
M = the total number of shiploads of cargo 
(movement requirements). 
NCOL = the required number of "good" schedules for 
ship j. 
NPICK = the minimum number of cargoes to be picked up 


by each good schedule. 
Vj = the dual price for ship j. 
uy = the dual price for movement requirement i. 
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step l 


Step 2 
Step 3 


Step 4 


step 5 


S21) 


Siz ) 


53) 


1 if movement requirement i is compatible 
with ship j 


O otherwise. 
For each i=1,...,M, set 
ug = uy*cyy + 99*(1 - cyy). 
Remove all movement requirements i with 
u,;>0 from further consideration, and let M' be 


the remaining movement requirements, 


i. e., those with u;sO. 


Sort uj; in an ascending order. 
Set KOUNT = O, P = {0}, and initialize the 
STACK. 
Put the movement requirements into the 
STACK in the same order as the sorted dual 
variables uj. 
While the STACK is not empty, do the 
following: 
Remove the movement requirement, say, k, from 
the top of the STACK. 
Calculate the completion time for path 
Bei stk js 
If this completion time s g, then do the 
following: 


a) Set P = PU {k}. 


1d 


b) Put the movement requirements not in the paen 
on top of the STACK in the order of decreasing 
Uy. 

cc) ‘Go to Step 5.2 

5.4) Otherwise, calculate the reduced cost for the 
schedule corresponding to path P. 

a) If the reduced cost is negative and the 
number of movement requirements in P r NPICK, 
then a good schedule is found. Set KOUNT = 
KOUNT + 1. 

If KOUNT s NCOL, remove node k from the STACK 
and go to Step 5.1. 
Otherwise, return to the master problem. 

b) If either the reduced cost is nonnegative or 
the number of movement requirements in P is 
less than NPICK, then remove node k from the 


STACK and go to Step 5.1. 


Note that Step 5.3.b dynamically generates the new 
level of nodes in the search tree. Moreover, the order in 
which the nodes are put on top of the STACK allows the 
search to progress in the depth-first manner with the 
movement requirement with the most negative cargo dual 
being searched first. In Step 5.4, the removal of node k 
from the STACK is equivalent to fathoming the node in the 


search tree. 
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IV. COMPUTATIONAL RESULTS 


The algorithm for generating schedules presented in 
Chapter 3 was implemented in VS FORTRAN on the IBM 3033AP 
computer at the W. R. Church Computing Center of the Naval 
Postgraduate School. The master problem is also implemented 
in the same manner and is fully described in the Master's 
thesis by Lima [Ref. 10]. 

For the computational experiments, we consider three 
problems which contain approximately the same number of 
lift assets and movement requirements as one would expect 
in a deployment of a small armed force such as that of 
Norway. Problem 1 has 48 movement requirements and 25 lift 
assets, Problem 2 has 48 movement requirements and 30 lift 
assets, and Problem 3 has 60 movement requirements and 30 
lift assets. The ports of embarkation (New York, Norfolk, 
Charleston, Jacksonville, Pensacola) and debarkation 
(Hamburg, Wilhelmshaven, Rotterdam, Antwerpen, Cherebourg) 
for the movement requirements are given in Table 4-1. Table 
4-2 and 4-3 give the distances between the various ports 
and the initial position of each asset. The speed for the 
assets are between 13 and 22 knots, and the probability 
that the ship is compatible with a given movement 


requirement is .75. 


IY, 


TABLE 4-1 : NUMBER OF MOVEMENT REQUIREMENTS 
BETWEEN POE'S AND POD'S. 
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TABLE 4-2 : DISTANCE BETWEEN POE 'SmalibaLo bis 
IN NAUTICAL MILES. 
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A. A HEURISTIC FOR SELECTING MOVEMENT REQUIREMENTS 

In Step 1, the algorithm discards those movement 
reguirements which have a positive dual value because of 
incompatibility with the ship under consideration. The 
remaining movement requirements are then selected for 
branching based solely on the magnitude of the associated 
cargo duals. Intuitively, this seems inefficient. Assume 
that the algorithm is constructing a new schedule for ship 
j}. Then those movement requirements which are in the 
schedule for ship j and those which are still unassigned at 
the end of the preceding master problem should be selected 
for branching first. Otherwise, movement requirements which 
have been assigned to other ships may be assigned to the 
new schedule for ship j also. This would create a "double 
coverage" for the same movement requirement. To create the 
least amount of double coverage a negative number is added 
to the cargo duals for the movement requirements which are 
already assigned to ship j and to those which have yet to 
be assigned to any ship. This negative number should be 
large enough to insure that when the movement requirements 
are ordered in Step 2 and put into the STACK in Step 3 and 
Step 5.3.b, the desired movement requirements are always 
above the others in the STACK. 

To implement the above heuristic, Step 1 of the 


algorithm in Chapter 3 is modified as follows: 
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Step 1 


a) For each 1 = 1,..., M set 
99 if Ciy = 0~ 
Uj if cj; = 0 and cargo i has been 


assigned to another ship l, 
where 1 = j, by the previous 
uj = < master problem iteration. 
= if cj3 = 0 and cargo i is 
unassigned or has been assigned 


to ship j by the previous 
master problem iteration. 


where 6 is a sufficiently large positive number. 
b) Remove all movement requirements 1 with 

u;>0 from further consideration, and let M' 

be the remaining movement requirements, i.e., 

those with uj;<0. 

The algorithm in Chapter 3 has been implemented both 
with and without the above heuristic for selecting movement 
requirements. Figure 4-4 shows that, based on the average 
CPU time to solve the three sample problems with various 
values for Tt, the algorithm with the heuristic clearly 
outperforms the one without. In fact, the algorithm with 
the heuristic is on the average 47% faster. Henceforth, the 


subproblem algorithm is implemented with the heuristic for 


selecting movement requirements. 
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Figure 4-4 


Comparison of total CPU time with and without 
heuristic implemented. 
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B. THE NUMBER OF SCHEDULES TO GENERATE FOR THE MASTER 
PROBLEM 


In general, it is unclear what is the optimal number of 
columns to generate during each cycle of the decomposition 
process. More columns mean more information for the master 
problem during each cycle. At the same time, more columns 
in the master problem also mean more CPU time for pricing 
out and other housekeeping operations. However, Figure 4-5 
demonstrates that it is advantageous to generate between 10 
and 20 schedules (columns) for Problem 2 when 7T = 29 days 
(the optimal duration). 

Table 4-6 provides an explanation for the above 
finding. The table lists the average number of nodes in 
the search tree for the two strategies: one column and ten 
columns. Although the average number of nodes in the search 
tree for the ten column strategy is more than the one 
column strategy in the early iterations, the information 
provided by the additional nine columns actually helps the 
master problem to find a solution quicker. Moreover, as the 
number of iterations grow, the one column strategy 
occasionally generates unacceptably large search trees in 
order to find one negative reduced cost column which could 
have been generated earlier if the ten columns strategy was 


employed. 


Zo 


CPU SECONDS IN SUBPROBLEM 


[2 


80 


48 MOVEMENT REQUIREMENTS 


SOR Sries 
OPTIMAL DEPLOYMENT DURATION 





10 20 50 
NUMBER OF SCHEDULES GENERATED 


Figure 4-5 


CPU time for subproblem versus 
number of schedules generated. 
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TABLE 4-6: AVERAGE NUMBER OF NODES GENERATED BY THE 
ALGORITHM FOR PROBLEM 2. 


One column per cycle Ten columns per cycle 
Iter Average # Iter Average # 
nodes 


nodes 

—. 

2S a 
pase | ae | ae |e 
po | we face |e 
paso [we [aes [os 
rao |e [| 
fame [| ue | | 
a 
SS 
paso | os [| 
pene| we | | 
re A 
awe| ae | | 
Pue| 7s | | 
aac [ was | | 
Psaco[ on | | 
pane] {| | 
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Figure 4-7 and 4-8 depict the trade-off between time 
spent in the master and subproblem for Problem 2. Figure 4- 
7 shows that by generating more columns per iteration one 
uses less CPU time for the subproblem. However, more 
columns per iteration also means more columns for the 
master problem which in turn leads to more CPU time for 
pricing out and other-housekeeping operations. In terms of 
total CPU time for the overall problem, Figure 4-8 shows 
that the 10 column strategy 1s better than the one column 
strategy when fr is less than or equal to the optimal 
duration (29 days). However, when T 1S bigger than 29 days, 
the one column strategy 1s superior. This is due partly to 
the fact that when +r is larger than the optimal duration, 
the problem 1S more relaxed, l1.e., it has a larger feasible 
region, and accurate dual information is not as critical as 
in the case when 7 is Smaller than the optimal duration. 

Table 4-9 summarizes the computational results for all 
three problems. It is clear from this table that the ten 
column strategy dominates the one column strategy for all 
three problems. However, it would be premature to make any 
conclusive recommendations based on the results presented 


here. 
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TABLE 4-9 : COMPUTATIONAL RESULTS FOR THE THREE PROBLEMS. 


PROBLEM 1 2 PROBLEM 2 a PROBLEM 3 


CPU CPU CPU CPU 
time time time time 


won| +f af ape 
Paverese [x | | ee | co | coe | pee 
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V. CONCLUSIONS AND FUTURE RESEARCH 


An algorithm for generating schedules for the crisis 
deployment problem is described. This algorithm is similar 
to the branch and bound algorithm for integer programs or 
the depth-first search technique for traversing directed 
graphs. The computational results in Chapter 4 demonstrate 
that on medium size problems the algorithm, when integrated 
with the master problem algorithm by Lima [Ref. 10], is 
effective in solving the subproblem of the crisis 
deployment problem. However, a more extensive computational 
study is required before any conclusive recommendations can 
be made. 

In addition to further computational studies, other 
possible areas for future research are listed below. 


(1) Allow movement requirements to arrive at POE's within 
a time window. 


(2) Allow movement requirements in partial shiploads. 
This would imply that a ship can pick up cargoes at 
two or more locations before delivering them to their 
discharging ports. 


(3) Investigate other branching strategies for the 
algorithm. For example, one can branch to the closest 
movement requirement instead of the movement 
requirement with the most negative dual value. 


(4) Investigate other combinatorial algorithms for the 
subproblem. 


(5) Consider possible delays due to random events such as 
attacks on ships or ports. 
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APPENDIX 
FORTRAN CODE FOR SUBPROBLENM 


SUBROUTINE SUBPR(U, XCOL, MDAY, M, N, NLOA, NPOE, NPOD, NSH, MR, TR, IT, A 
MEGHP AT, 1H90e) SHIP, KTIME, Ki SPD) SEQ, LSEO, CTSHIP, GEINODE, NPICK, NCOL) 


e THIS SUBROUTINE APPENDS NEW (GOOD) COLUMNS TO THE A MATRIX IN THE ° 
° MASTERPROBLEN. © 
e e 
® KEY VARIABLES COMING FROM THE MASTER PROBLEM: ° 
e e@ 
zee Us DUAL VARIABLES. ° 
e - MDAYs MAXIMUM DURATION OF ANY SCHEDULE. ° 
= = t A MATRIX. ° 
© - COMPAT: SHIP COMPATIBILITY WITH SHIPLOADS OF CARGO. ° 
° - IB: INDEX SET FOR THE BASIS. ° 
° - XB: COLUMNS IN THE BASIS ° 
2 > NPICK: HINIMUM NUMBER OF SHIPLOADS OF CARGO TO BE PICKED BY ° 
© = NCOL: NUMBER OF SCHEDULES TO GENERATE FOR EVERY CALL TO THE ° 
, SUBPROBLEN. ° 
e ® 
® KEY VYARIBLES USED IN THE SUBPROBLEM: ° 
e ° 
°e - VIND: INDEX SET OF re SORTED DUAL VARIABLES. e 
=) - LOAD: SHIPLOAD NUMBE : 
zoe PATH: SEQUENCE OF SHI PLOAD NUMBERS FOR THE CURRENT SCHEDULE ° 
> = STACK? SHIPLOAD NUMBERS IN THE STACK. ° 
wae CURLED: CURRENT SHIPLOAD NUMBER. ° 
woe PRED PREDECESSOR OF SHIPLOAD NUMB . 
fae TLE COMPLETION TINE FOR GENERATED T SCHEDULES. ° 
© - LENGTH: LENGT OF CURRENT PATH. id 
® e 
e e 
@e@rerwrwee mi ew ew Hwee ewe ewe wee ewe Hew ew ew ewww wwe wee ew ewewwewee ewe ewer e2eweewewewewee2xe w@ w® e@ we ew e@ ww ew t 


IMPLICIT DOUBLE PRECISION (A- =H, O- 25) TPEGER  r- =H) 
PARAMETER( MM #100, NN = 200 = JJ 000 
DIMENSION KEOL CHM, UCHHD, OUCH, VCRND, XBUMAD, SACHA, SPDCHM), TBC MM) 

A MM 
INTEGER VIND( MM), PRED(O1IJ), LOAD(O1:JJ), TIME(O:JJ), STACK(O1JJ), TOP, 

& MR(100, KK), TR(15, 15), IT( 30, 5), FROLD, TOLD, PATH(0:4M), CURLD, COUNT, — 

& LENGTH, LASTND, SHIP. TT, MLNGTH, KTIME(NN), CTSHIP, CTBACK(OHHD, 

& CTSON(O1MM), SEQ(NN, MM), LSEQ(NN), CTNODE’ 


COMMON /UNITS/ NIN, NOUT 
C INITIALIZE 
CTNODE s 


0 
LIMIT = N 
CTSHIP # O 


360 CONT 


NNEG © NLOA 
MINRC # O.ODO 


eeeeecese eweewewereeeewe ew ewww ew egw ei ew weeeeeeweew wv’ vw ew@eeeeewewreeewrevrewerewewwwownrerrr2r 2 2°2 222727 °°" 


© SORT THE DUAL VARIABLES : 
Sm wen ween een eww wen ene wn eH eee ee ew He Hew em ew He ew ew KH He ee ee eH em ee em em er em emer err ewrroererrrrrrrrrr: 
DO 8&8 I # 1,NLOA 
8 UCT) = UCT) 
GC SREURISTIC 
C ASSIGN LARGE NEGATIVE DUAL VALUE TO SHIPLOADS NOT ALREADY PICKED 
© BY GIHER SHIPS 
TEGenN. GE. (2eNLOA * NSH © 1) ) THEN 
DOy2ZO1 1. * 1,8 
DO 202 J 2#= 1,4 
IF( (IB(J) .GT. (2*NLOA * NSH)) . AND 
& (ABC I) .GT. DO) ) THEN 
SA(I) # SACI) + A(T, IB(J)) 
END IF 
202 CONTINUE 
IF( SACI) .LT. 1.0D0O ) UU(I) #® -2.0D0 
201 CONTINUE 
END IF 
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C REMOVE SHIPLOADS NOT COMPATIBLE WITH SHIP FROM FURTHER CONSIDERATION. 


DO 9 I = 1,NLOA 
9 IF( .NOT. COMPAT(I,SHIP) ) UU(I) = 99,0D0 
DO 20 I = 1, NLOA 
MIN = 0.1D-6 
COUNT = O 
IND « 
DO 30 J = 1,NLOA 
UUcy) .LT. MIN ) THEN 
MIN = UU(J) 
IND = J 
COUNT = 1 
END IF 
30 CONTINUE 
IF( COUNT .EQ. O ) THEN 
NNEG = I - } 
GO TO 40 
END IF 
V(I) = MIN 
VIND(I) = IND 
UUCIND) = 99.0D0 
20 CONTINUE 
40 CONTINUE 
@re ee eke eee eee ee me me ee em me we ew Me me wm mem we me wm wm mm we ewe wm wm wee wm we we ewe www ew ew ewe ew eee eee e 
* DEPTH FIRST SEARCH 
re eee ee ee ee ee ) 
GC INTTIALIZE 
DO 350 I #0,JJ - 1 
TIME(I) = O 
STACKUI) <i @ 
LOAD(I) = O 
PRED(I) = O 
350 CONTINUE 
DO 380 I = 1,4MM 
XCOL(1I) = 0.DO 
PATH(I-1) = O 
CTBACKCI-1) «= 0 
CISGn (tel) «0 
380 CONTINUE 
LENGTH = O 
LASTND = } 
TOP = O 
CURLD «= 0 
FROLD = O 
TOLD = O 
C CREATE ALL NODES OUT OF THE SOURCE AND PUT THEM IN STACK 
DO SO I «= NNEG,12,-12 
LOAD(NNEG - I * 2) «= VIND(I) 
PREO(NNEG - £ +02) end 
STACK(NHEG. = £ *« 1) = tNNEG = 1o02) 
TOP = TOP 


gee 
LASTND #® LASTND + 1 
SO CONTINUE 


C MAIN LOOP TO SEARCH FOR FEASIBLE SCHEDULES 
100 CURLD = STACK(TOP) 
C SPECIAL CASE WHEN ONE SHIP CAN PICK ALL SHIPLOADS OF CARGO 


IF( LENGTH .E@. NNEG ) THEN 
RCOST = QO. ODO 
DO 91 I = 1,LENGTH 
91 XCOL(LOAD(PATH(I))) =» 1.0D0 
RCOECNLOA * SHIP) #= 1.0D0 
DOL Tiel 1,0 
iii IF ( XCOL(I) ~-EQ,. 1.0D0 RCOST s aoe * UCT) 


IF CRCOST .GT. <1. 0D-47 RCOST =» O. 
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aaa (RCOST LT. 0.O0DO) ) THEN 
DO 389 I = i,m 
389 ACI,N) © XCOL(I) 
DO 122 Jei, LENGTH 
22 SEQ(N, J) © LOAD(PATH(J)) 
LSEQ(N) #» LENGTH 
RETURN 
END IF 
END IF 


C IF STACK EMPTY RETURN TO MASTERPROBLEN IF GOOD SCHEDULE FOUND. 


IF( CURLD .EQ,. O ) THEN 
IF( MINRC .LT. 0O.ODO ) RETURN 


C IF NO GOOD SCHEDULES FOUND RUN SUBPROBLEM AGAIN FOR ANOTHER SHIP. 


CISHIP = CTSHIP * 1 
SHIP SHIP *) 2% 


IPCSHIP .EQ. NSHos 1) SHIP = ft 
tF¢ CTSHIP . NE. Nsw 2? GO TO sit 
RETURN 

END IF 


TEC PREDCUCURLD) .~EG. 1) THEN 
LASTND = LASTND - LENGTH 


LENGTH « O 
END IF 
IF( CTSON(LENGTH) .E@. CTBACK(LENGTH) ) THEN 
DO 51 I «© 1, LENGTH 
IF( PRED(CURLD) .EQ. PATH(I) ) THEN 
LASTND = LASTND - LENGTH + I 
LENGTH = I 
GO TO 539 
END IF 
51 CONTINUE 
END IF 


339 PATH(LENGTH + 1) = CURLD 


FROLD © LOAD(PRED(CURLD) ) 
Ee . LOAD(CURLD) 


C CALCULATE COMPLETION TIME FOR CURRENT SCHEDULE. 


CALL TIMEFC(NPOE, NPOD, MR, TR, IT, FROLD, TOLD, SHIP, NSH, NLOA, TT, SPD) 
TIME(CURLD) © TIME(PRED(CURLD)) * TT 


C IF FEASIBLE INCLUDE SHIPLOAD NUMBER IN PATH, 
TEC TIMECCURLD) . CE. MDAY ) THEN 


CTBACK(LENGTH) # CTBACK(LENGTH) + 1 
LENGTH # LENGTH + 1 
CTSON(LENGTH) # O 


C PUT INTO THE STACK ALL SHIPLOAD NUMBERS NOT ALREADY IN THE PATH. 
DO 60 1 = NNEG,1,-1 


DO 70 J = 1, LENGTH 
70 TP ClVINDCI) 2 EG. LOAD (PATH CIO GO To 60 
LASTND « LASTND «+ 1 
LOAD(LASTND) ® VIND(I) 
PRED(LASTND) = CURLD 
TOP « TOP + 
CTSON(LENGTH) #® CTSON(LENGTH) + 21 


60 CONTINUE 
DO 8O I _= LASTND, (LASTND - CTSON(LENGTH) + 1),-1 
STACK(TOP) = I 
TOP © TOP - 
80 CONTINUE 
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¢ 


TOP 


TOP 


ELSE 


LASTND #® LASTND - 1 


RCOST # QO. ODO 


CTSON (LENGTH) 


C CHECK IF GOOD SCHEDULE FOUND. 


90 


110 


388 


222 


400 


130 


C 


DO 90 I #®.1,LENGTH 
XCOL(LOAD(PATH(I))) 


XCOL(NLOA * SHIP) 


DO 110 J -2 i,n 
IF( XCOL(I) .EQ. 


IF(RCOST .GT. 
IF( RCOST .LT. 


0. ODO 


IF ( CTBACK(LEN 


N = N + 
DO 388 I 
A(I,N) ® 


XTIME(N) 


DO) 542. ° 
SEQ(N, J) = 


LSEQ(N) L, 


IF(¢ RCOST . 
MINRC s 
K = N 
LE 


XC 


i, 


END 
ENDGLe 


DO 400 I = 1,LE 
KCOL(LOAD(P 
XCOL(NLOA ¢ 

CONTINUE 


ELSE 


DO 130 . 

KCOLY LOAD 

XCOL(NLOA 
END IF 


1.0D0 ) 
~1.0D-4) 


LEN 
PA ATHCL 
* (SHIP 


1. ODO 
1. ODO 


RCOST = RCOST + U(T) 
RCOST # O.ODO 


LENGTH .GT. NPICK ) 


» AND. 


GTH) Eee G, THEN 
M 

oL«1) 

TIME(PRED(CURLD) ) 


LENGTH 
LOAD(PATH(J)) 


ENGTH 


LT. MINRC ) 
RCOST 


THEN 


NGTH 
ATH(I))) 
SHIP) » 


~~ G) 
~~ 
a~r 


IF NO GOOD SCHEDULE FOUND BACKTRACK. 


END 
TE < 
CTB 
GO 


END 


eh ae dala 

GO TO 10 

If 

N .GE. LIMIT*NCOL) 
ACK(LENGTH) = O 
TO 100 


» CTBACK(LENGTH) 


td 


1 


RETURN 
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THEN 


Sr a ed TIMEFC(NPOE, NPOD, MR, TR, IT, FROLD, TOLD, SHIP, NSH, NLOA, TT, 


—see-eoweeeweweweree2ee2e2e2@ e222 @@@eww @ wT ww w@BFw=eawweweeeeewweewreeweeseeeeeeeweqee ewes eeeaoeqeeqeoe® @ = 


* THIS SUBROUTINE CALCULATES THE TIME REQUIRED TO PICK UP AND DELIVER 
® A SHIPLOAD OF CARGO IN WHOLE DAYS. 


IMPLICIT DOUBLE PRECISION ( A-H,O-Z ), INTEGER ( I-N ) 
PARANETER( KK = 2, MM #100, NN #® 2000 ) 


DIMENSION SPD(MM) 
INTEGER MR( 100, KK), TR(15, 15), IT(30, 5S), TT, TOLD, SHIP, FROLD 


TT = O 
C CALCULATING THE TRAVEL TIME 
IF(FROLD .EQ. OQ) THEN 


TT = IDNINT(C(CIT(SHIP, MR(TOLD, 2)) * TR(MR(TOLD, 1), MR(TOLD, 2)))/ 
(24. © SPD(SHIP))) 


& 
ELSE 
TT = IDNINTCCTR(MRC(TOLD, 1), MRC FROLD, 2)) + TRC(MRC(TOLD, 1), 
& MRCTOLD, 2)))/(24. © SPDC(SHIP))) 
END IF 
RETURN 
END 


Sy 


oe 


9. 


10. 
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